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Chapter 1 
Introduction 


This is an semi-annual report that summarizes the work done on Grant 
No. NAG2-542 during the period 1 September 1988 to 28 February 1989. 
There are four areas of study to be summarized: 

1. Bistatic scattering measurement studies for a compact range. This 
work is being conducted by T. H. Lee and W. D. Burnside under Ohio 
State University project number 721432. 

2. Target signature modeling for test and evaluation hardware in the 
loop situations. This work is being conducted by R. J. Marhefka 
under Ohio State University project number 721446. 

3. Aircraft and missile code modification studies. This work is being 
conducted by R. Rojas under Ohio State University project number 
721447. 

4. SATCOM antenna studies on aircraft. This work is being conducted 
by R. J. Marhefka and D. Bensman under Ohio State University 
project number 721711. 

An overview of the research accomplished is described in the separate 
chapters below. 
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Chapter 2 

Bistatic Measurement Study 


This is a brief review of the portion of the grant pertaining to bistatic 
measurement study over the period September 1, 1988 to February 
28, 1989. The summary described here is to be presented at the 
AMTA meeting in August entitled “Imaging Processing of Bistatic 
Scattered Fields Obtained in a Compact Range,” by T. H. Lee and 
W. D. Burnside. 


2.1 Introduction 

It is very important in diagnostic radar applications that one can 
accurately distinguish various dominant scattering terms associated 
with specific features on a target based on its scattering characteris- 
tics. Normally, the backscattered fields of the target are collected and 
processed to obtain the down-range and cross-range scattering cen- 
ters [1], The backscattered fields can be collected either in a far-field 
range or in a compact range. On the other hand, bistatic scattering 
of a target has been getting a lot of attention recently. Consequently, 
the determination of the dominant scattering centers using bistatic 
scattered fields of a target is also very important. 

In this report, a general image processing technique which utilizes 
the bistatic scattered fields collected in a compact range is presented. 
The transmitting and receiving antennas can be either two compact 
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ranges or one compact range and a horn antenna. The compact range 
reflector can be either focussed or defocussed so that a near field sit- 
uation can be simulated. The bistatic scattered fields are collected 
as a function of frequency and the angle of rotation of the target. 
Then, cross-range and down-range images of the target are deter- 
mined. Experimental results are presented to validate this image 
processing technique. 


2.2 Approach 


Conventional image processing of backscattered fields for a target 
uses efficient algorithms such as a Fast Fourier Transform to deter- 
mine the target image. This requires that the antenna be located a 
large distance from the target [2] and that the target be rotated only 
through small angles. However, these assumptions are not necessary 
in general. For example, it has been shown that one can use defo- 
cussed compact range system to simulate a near field illumination of 
the target [3]. In order to determine the image of the target from 
the collected backscattered fields, a more general image processing 
approach has to be implemented [2,3]. The key factor in general im- 
age processing is the knowledge of the phase history of the scattering 
center associated with a given illumination. The backscattered fields, 
which are normally collected as functions of frequency and target an- 
gular rotation, are processed by multiplying the received response by 
a phase factor which is the complex conjugate of the phase associated 
with the path length from the source of illumination to the potential 
scattering center and then to the receiver. This data is then summed 
to determine the image response. Using this approach, the response 
at a true scattering center is more enhanced than those which are not 
true scattering centers. 

The same processing technique can be used to determine the scatter- 
ing centers of a target from its bistatic scattered fields. Consider the 
bistatic scattering measurement configuration shown in Figures 2.1 
and 2.2. The system can be either two compact ranges as shown in 
Figure 2.1 or one compact range and one horn antenna (or other an- 
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Figure 2 . 1 : Top view of bistatic scattering measurement in a compact range. 
(Transmit and receive compact ranges.) 

tennas) as shown in Figure 2 . 2 . For a target located in the test zone, 
the bistatic scattered fields, E,(f m , 9 n ), are collected for a series of 
frequencies, / m , and a series of angles of rotation, 0 n . The image of 
the target, /( 7 *, 7 Z ), is determined by 

/(7..7.) = "if "if (2.1) 

m—1 n=l 

where 

*( 7 ., 7 .,/rn,*„) = — (-1 + s 2 ). ( 2 . 2 ) 

c 

In Equations 2.1 and 2 . 2 , 8 \ and s 2 are the path lengths from the 
transmitting antenna to the point ( 7 ,, 7 Z ) on the image plane and 
from (7*, 7 *) to the receiving antenna, / m is the m’th frequency and 9 n 
is the n’th angle of rotation. It is apparent that if the point, (7*, 7*), 
is a scattering center, the bistatic scattered field, E t , should have a 
phase factor, e~* 2 *, which is the complex conjugate of the phase factor 
shown in Equation 2 . 1 . Consequently, the resulting image, /( 7#, 7*), 
will peak up at the image location. On the other hand, if (7,., 7*) 
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Figure 2.2: Top \ iew of bistatic scattering measurement in a compact range. 
(Transmit compact range and receive horn antenna.) 

is not a scattering center, the image will not be as significant. For 
each point on the image plane, one can determine the path lengths, 
3 i and s 2 , and use Equation 2.1 to determine the complete image of 
the target. The frequency behaviour of the bistatic scattered fields 
provides the down range information; whereas, the angle responses 
provide the cross range location of the scattering centers. 

Since the image processing method given by Equation 2.1 only re- 
quires the knowledge of the path lengths, Sl and s 2 , the compact 
range reflector can be a defocussed system so that one can simulate a 
near field situation. However, one has to determine the ray trajectory 
in order to find the point of reflection on the compact range reflector. 


2.3 Results 


A set of bistatic measurement has been performed at The Ohio State 
University ElectroScience Laboratory compact range facility. The 
configuration of the measurement is shown in Figure 2.3. Spheres 
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Figure 2.3: Bistatic scattering measurement of spheres in OSU compact 
range. 
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with diameters of 6 ”, 3.25” and 1 ” were placed on a styrofoam plat- 
form which was then positioned at the center of the target zone. The 
compact range reflector has a focal distance of 12 ’ and the center of 
the target zone is at 24’ from the vertex of the reflector with a height 
of 60”. The second antenna for the bistatic scattering measurement 
was located 193” from the center of the target zone in the direction 
perpendicular to the axis of the compact range reflector. The height 
of this antenna was the same as the center of the target zone. Several 
combinations of spheres were used as the target for these measure- 
ments. The target was rotated from -30° to 30° with respect to the 
axis of the compact range reflector, and bistatic scattered fields were 
collected uniformly across the 2 to 18 GHz frequency band. Since 
the compact range has a capability of simultaneously collecting data 
from 2 different channels, the backscattered fields were also collected 
at the same time. Two image processing results based on these mea- 
surements are presented in this section. 

The first combination of the spheres used as the target is shown in 
Figure 2.4. Note that a 6 ” sphere is located at X and a 3.25” sphere 
at N. The image of this target obtained by processing the bistatic 
scattered fields is shown in Figures 2.5 and 2.6. The backscattered 
fields can also be processed by using Equation 2.1 by simply assuming 
that the two compact range reflectors shown in Figure 2.1 are co- 
located. The resulting image is shown in Figures 2.7 and 2.8. From 
Figure 2.4, the center of the 6 ” sphere is at 7 X =8.485” and 7 I = 0 . 0 ” 
and the center of the 3.25” sphere is at =-4.243” and 7 * =-12. 728”. 
However, the scattering centers will not be at these locations because 
one has to take the radius of the spheres into account. The images 
shown in Figures 2.5 to 2.8 clearly indicate the scattering centers of 
this target. One also notices the movement of the sphere scattering 
centers from these images. This is because the target was rotated 
with respect to the center of the target zone which results in the 
movement of the scattering centers in that the sphere reflection point 
moves with the angles of incidence. In addition to these scattering 
centers, there are other peak responses that appear in these images. 
These are secondary scattering centers caused by the creeping waves 
associated with each sphere. 
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Figure 2.4: Positions of the 6” and 3.25” spheres for the bistatic scattering 
measurement. 
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Figure 2.5: Surface plot of the images for the 6” and 3.25” spheres obtained 
from the bistatic scattered fields. 
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Figure 2.7: Surface plot of the images for the 6” and 3.25” spheres obtained 
from the backscattered fields. 
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Figure 2.9: Positions of the 6”, 3.25” and 1” spheres for the bistatic scat- 
tering measurement. 


The second combination of the spheres used as the target is shown in 
Figure 2.9. Note that a 6” sphere is located at V, a 3.25” sphere at J, 
and a 1” sphere at L. The resulting images obtained from the bistatic 
scattered fields and backscattered fields are shown in Figures 2.10 
to 2.13. Again, these images clearly show the scattering centers of the 
target and also the secondary scattering centers due to the creeping 
waves propagating around the spheres. 
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Figure 2.10: Surface plot of the images for the 6”, 3.25” and 1” spheres 
obtained from the bistatic scattered fields. 
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Figure 2.11: Contour plot of the images for the 6”, 3.25” and 1” sphen 
obtained from the bistatic scattered fields. 
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Figure 2.12: Surface plot of the images for the 6”, 3.25” and 1” spheres 
obtained from the backscattered fields. 
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2.4 Summary and Conclusions 


An image processing technique using the bistatic scattered fields ob- 
tained in a compact range has been presented. Two compact ranges 
or one compact range and a horn antenna can be used to perform 
the bistatic scattering measurements. The antennas can be either fo- 
cussed or defocussed. The bistatic scattered fields are then used to 
find the image, and consequently, the scattering centers of the target. 
The same processing approach can also be U6ed to find the image of 
the target from backscattered fields. The processing technique does 
not require the target to be at far field of the antennas. Experimental 
results have been used to determine the images of various targets via 
this processing technique. Very good results have been obtained and 
presented. 


2.5 Acknowledgement 

Numerical results presented in this report were made possible by a 
Grant (Project 355) from The Ohio Supercomputer Center, Colum- 
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Chapter 3 

Target Signature Study 


This is a brief review of the portion of the grant pertaining to target 
signature modeling over the period September 1, 1988 to February 
28, 1989. 

3.1 Introduction 

The objective of this part is to determine and implement an effec- 
tive method to integrate the use of UTD electromagnetic prediction 
codes for simulation of radar cross section into test and evaluation 
situations. UTD scattering code predictions can replace costly mea- 
surements. The signatures are used in the test loops to evaluate how 
effectively radar can follow targets. In this situation, it is necessary 
that the data be efficiently obtainable. 


3.2 Scattering Code Overview 


There is much scattering code development going on in the aircraft 
industry. These codes are primarily physical optics codes based on 
tiled or patched surfaces. They are used for component analysis and 
require large amount of computation times to do the integrations. 
Because they are based on physical optics, the specular points are 
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not explicitly defined. In addition, the some of the codes can be 
proprietary and/or very expensive. 

A physical optics code is being developed by the U. S. Army Missile 
Intelligence Agency and Northrop called MISCAT. It is available in- 
expensively, however, it has the characteristic mentioned above and 
in addition has a more basic shape input capability. A new code is 
being developed called PGSCAT, which is a patch code for more gen- 
eral model descriptions, but again, it is PO based. Syracuse Research 
Corp. has a code that over comes the PO integration by using a scat- 
tering center description. There are no higher order terms involved 
and it uses basic shape models. 

A different class of scattering codes are centered around the BRL 
modeling code. The BRL code is a CSG modeller for graphical rep- 
resentation of shapes. ERIM’s SRIM code uses the BRL brute force 
ray shooting methods to reproduce synthetic aperture radar (SAR) 
images. It provides approximations to the fields as projections of the 
rays onto pixel planes. The ray shooting technique does find dihedral 
and trihedral speculars but it is extremely time consuming. GTRI 
has a RCS code similarly centered around the SAR concept. Univer- 
sity of Illinois has taken the BRL ray tracing code and added a field 
capability to it. The resulting fields lack edge effects and is slow since 
it is based on PO and brute force ray tracing. 

Another class of scattering codes deal with the near zone fusing prob- 
lem. NWC has integrated the Ohio State University’s NEC-BSC code 
into its applications for this purpose. 

The suggested approach for test and evaluation purposes, is the use 
of a UTD based code such as the RCS - Basic Scattering Code [4]. 
It explicitly defines the specular, edge, and vertex diffraction points, 
as well as some higher order terms. It uses efficient ray tracing tech- 
niques based on the chosen geometry models. It does not need inte- 
grations since it uses closed form diffraction coefficients. This allows 
the flexibility to store the characteristics of the scattering centers of 
a target and then add the field information versus angle or frequency 
as they are needed. Since test and evaluation hardware in the loop 
situations need efficiency and flexibility of operation, this type of ap- 
proach will be beneficial for large sized targets and the many angles of 
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approach needed. In addition, it has a bistatic scattering capability 
that might be useful in the future. 

A study of the various UTD mechanisms that will be most important 
for modeling in the test and evaluation situation has been under- 
taken. Their theoretical characteristics have been studied from the 
point of view of obtaining the best way to determine, calculate and 
store each term individually. The UTD radar cross section prediction 
code, RCS-BSC version 2, has been modified so that its write algo- 
rithms for the UTD mechanisms calculated will better accommodate 
a volumetric output of the individual fields. 

The RCS-BSC version 2 corner diffracted field section as been modi- 
fied to allow access of an individual corner diffracted field instead of 
an edge in order to study which is the best way to store the fields 
in the data base. Other field sections will be modified as they are 
integrated into the data base. 

An algorithm to access and plot the three dimensional field infor- 
mation in the data base has been started. This will allow a visual 
presentation of the fields in order to help determine the minimum 
sample rate for the different individual types of UTD mechanisms. 


3.3 Storage Requirements 

In order to better understand the size requirements of the data base 
needed for test and evaluation purposes the following simple examples 
illustrate the various possible configurations. It is assumed that the 
this information is for backscatter at a single frequency at this point. 
It is also assumed that a complete volumetric view of the object under 
consideration is needed with all polarization information. 

For a full up unprocessed scattered field data base: 

• Opts ~ 1801, <j>pts = 3601, no. pols . == 4 

• Complex number = 8 bytes, step size = 0.1° 

• Storage required = 207 megabytes 
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For a full up processed data base, that is with a median average 
representation of the fields: 

• Opts = 91, <f>pts = 181, no. pols . = 4 

• Complex number = 8 bytes, step size = 2° 

• Storage required = 0.527 megabytes 

The median average of the fields is not necessarily desirable, however, 
for test and evaluation purposes since the scintillation of the fields is 
useful to test radar performance. In addition, the processed data base 
can not be used for imaging since the phase information is lost. 

For individual dominant terms only, it can be assumed that the sam- 
pling rate can be greatly reduce so that: 

• Opts = 37, 4>pts = 73, no. pols . = 4 

• Complex number = 8 bytes, step size = 5° 

• Storage required per term = 0.086 megabytes 

This implies that approximately 10-20 terms can be used and still 
have the storage and more field information content than the pro- 
cessed data base. A few thousand terms can be stored and have the 
same information as the full up unprocessed data base. Actually some 
of the storage would be needed for several frequencies. This, however, 
appears to be the best alternative. 

One more possible solution has been investigated. This entails stor- 
ing scattering centers only and then post processing the fields. The 
storage of the scattering centers locations only turns out to be very 
complicated since the terms can be shadowed or they can move ver- 
sus incident and observation direction. In addition, there will be a 
trade off between storage and speed of retrieval which will reduce its 
desirability for real time test and evaluation purposes. 

In summary: 

• Full up unprocessed data base requires large amounts of storage. 

• Processed data base loses local variations of signals. 

• Processed data base can not be used for imaging. 
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• Individual dominant terms can be reconstructed to give the same 
information as full up unprocessed data base with less storage. 

• Need for multiple frequency information further compounds full 
up data storage, 

• Individual terms data base will still greatly reduce storage re- 
quirements. 

• One, two or three dimensional imaging possible. 


3.4 Further Work 

Test cases, that will represent meaningful situations for test and evalu- 
ation purposes, are being determined to validate the algorithms being 
developed. Simple scattering shapes will be used such as plate and 
cone frustums to test the basic ideas. A model of a QF-100 aircraft 
has been suggested by K. Oh and C. L Yu at Pacific Missile Test 
Center, since they have measured data on a one fourth scale model. 
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Chapter 4 

Aircraft Code Study 


This is a brief review of the portion of the grant pertaining to mod- 
ifications to the aircraft and missile antenna code over the period 
September 1, 1988 to February 28, 1989. 


4.1 Introduction 

The purpose of this effort is to modify the aircraft and missile antenna 
code [5] to provide simple models for rectangular microstrip patches 
and to add a smoothing algorithm for radiation pattern presentation. 
In addition, a study to implement the aircraft code on a personnel 
computer will be initiated. 


4.2 Smoothing Algorithm 

The existing Ohio State University Aircraft code (NEWAIR3) [5] has 
been modified so that any spikes in the radiation patterns are removed 
by a smoothing algorithm. The input data is as follows: 

EX: 

LSMOOTH, NAVG, SEN 
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where 


LSMOOTH = T, spikes are removed from the output data 
F, the output data is not modified 


NAVG — $ of data points used to compute an average. 

It is always even, usually NAVG=4. 

SEN = # of standard deviations that is allowed before a 

spike is removed. SEN does not have to be an integer. 

4.3 Patch Antenna 


Another task that is near completion is the addition of a simple model 
for a rectangular microstrip patch. The patch is modeled by a pair 
of short magnetic monopoles which is accurate when the lowest order 
(1,0) mode is excited in the patch. The input data of the NEWAIR3 
code (in the SG: command) has been slightly modified to input the 
parameters for the rectangular patch. 


4.4 PC Study 


A study to implement the NEWAIR3 on a personnel computer has 
been initiated. Additional memory and software has been obtained 
to start the task on an IBM compatible clone. This work will be 
continuing during the next phase. 


25 



Chapter 5 

SATCOM Antenna Study 


This is a brief review of the portion of the grant pertaining to satellite 
communication antenna (SATCOM) modeling over the period from 
the start of this part December 1, 1988 to February 28, 1989. 


5.1 Introduction 

The purpose of this effort is to study the performance of SATCOM 
antennas on aircraft. The potential locations of the antenna system 
for best pattern and polarization coverage is being assessed on an P3C 
aircraft. A circular polarized antenna composed of crossed dipoles is 
being considered. The aircraft is being modeled using a computer 
code based on the Uniform Geometrical Theory of Diffraction (UTD) 
which is the NEC-BSC V3.1 [6]. 

5.2 Batwing Antenna 

The first phase has been to analyze the present candidate UHF SAT- 
COM antenna system and its location on a P3C aircraft using the 
NEC-BSC. Naval Air Test Center has provided the antenna system in- 
formation which is a Dome and Margolin DM 1501341 (Batwing) air- 
borne UHF satellite communication antenna as shown in Figure 5.1. 
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A simple crossed dipole has been used in the code with the corre- 
sponding dimensions given in the figure. All the results have been 
referenced to a circular polarized isotropic radiator. The calculated 
result in Figure 5.2 at 244 MHz can be compared with the measured 
result supplied by the manufacturer in Figure 5.3 for the antenna on 
a 8 foot ground plane in the principal elevation cut for the antenna. 
The results compare very well. In order to test the NEC-BSC’s rep- 
resentation of the antenna for a pattern taken along the length of 
one of the dipole arms, the pattern at 300 MHz in Figure 5.4 can be 
compared with a method of moments result [7] in Figure 5.5. Again 
quite good agreement is obtained. 


5.3 Model Validation 

Naval Air Test Center has supplied extensive model measurements 
independently conducted by Boeing on a 1/17 scale model and Lock- 
heed on a 1/10 scale model. The Boeing data is used to validate 
the model here. The antenna is located as illustrated in Figure 5.6, 
which also shows the computer model used to generate the results. 
The calculated results for 300 GHz are compared with measurements 
for the roll plane in Figures 5.7 and 5.8, for the azimuth plane in Fig- 
ures 5.9 and 5.10, and for the elevation plane in Figures 5.11 and 5.12 
all for right hand polarization. The cross polarized fields are com- 
pare for the roll plane in Figures 5.13 and 5.14, for the azimuth plane 
in Figures 5.15 and 5.16, and for the elevation plane in Figures 5.17 
and 5.18 all for left hand polarization. Notice that the roll plane 
results compare very well. The azimuth and elevation results, how- 
ever, have about a 10 dB difference near the horizon. This will be 
investigated further, but overall the agreement is good. 


5.4 Further Work 


In the next phase, a study of alternative antenna locations and possi- 
bly alternative systems to provide a minimum of 330° azimuth, from 
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OUTLINE DIMENSIONS 

inches (Centimeters) 






Figure 5.2: Calculated Batwing right hand circular polarized antenna pat- 
tern at 244 MHz on an 8’ ground plane using the NEC-BSC. 
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Figure 5.3: Measured Batwing right hand circular polarized antenna pat- 
tern at 244 MHz on an 8’ ground plane. (Pattern supplied by Dome and 
Margolin.) ^0 



Figure 5.4: Calculated Batwing right hand circular polarized antenna pat- 
tern at 300 MHz on an infinite ground plane using the NEC-BSC. 
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Figure 5.5: Calculated Batwing right hand circular polarized antenna pat 
tern at 300 MHz on an infinite ground plane using the ESP code. 
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Figure 5.6: Geometry of the model of the P3C aircraft code used in the 
NEC-BSC showing the location of the antenna. 
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Figure 5.7: UTD calculated roll plane pattern for batwing antenna on a 
P3C for right hand circular polarization in the roll plane at 300 MHz. 


34 



ORIGINAL PAGE IS 
OF POOR QUALITY 


■■g 

H9S 

mm 

ay 


U 


—O fn i a s bw. 

n»t mom 

tt- _UHFS ATCOM 2 
m> STA . 91 Q 

axiWr 


hzd 



Isotrooe ■ -10.7 dbic 


l*T S’ 




*< noo 

OOOTANT AMQLf ♦ - 27 Q 

loumoATKm « # < k. c # « l 

u 


RHC“ 


Figure 5.8: Boeing’s measured roll plane pattern for batwing antenna on a 
P3C for right hand circular polarization in the roll plane at 300 MHz. 
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Figure 5.9: UTD calculated azimuth plane pattern for batwing antenna on 
a P3C for right hand circular polarization in the roll plane at 300 MHz. 
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Figure 5.10: Boeing’s measured azimuth plane pattern for batwing antenna 
on a P3C for right hand circular polarization in the roll plane at 300 MHz. 
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Figure 5.11: UTD calculated elevation plane pattern for batwing antenna 
on a P3C for right hand circular polarization in the roll plane at 300 MHz. 
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Figure 5.12: Boeing’s measured elevation plane pattern for batwing antenna 
on a P3C for right hand circular polarization in the roll plane at 300 MHz. 
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Figure 5.13: UTD calculated roll plane pattern for batwing antenna on a 
P3C for left hand circular polarization in the roll plane at 300 MHz. 
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Figure 5.14: Boeing’6 measured roll plane pattern for batwing antenna on 
a P3C for left hand circular polarization in the roll plane at 300 MHz. 
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Figure 5.15: UTD calculated azimuth plane pattern for batwing antenna 
on a P3C for left hand circular polarization in the roll plane at 300 MHz. 
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Figure 5.16: Boeing’s measured azimuth plane pattern for batwing antenna 
on a P3C for left hand circular polarization in the roll plane at 300 MHz. 
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Figure 5.17: UTD calculated elevation plane pattern for batwing antenna 
on a P3C for left hand circular polarization in the roll plane at 300 MHz. 
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Figure 5.18: Boeing’s measured elevation plane pattern for batwing antenna 
on a P3C for left hand circular polarization in the roll plane at 300 MHz. 
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zenith to 10° above aircraft horizon will be undertaken. The antenna 
system must provide performance to maintain the satellite link in ar- 
eas of high signal fading to multipath effects. This will primarily be 
accomplished with the use of the same antenna mentioned above by 
exploring alternative locations around the aircraft. 
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